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Perturbative solution to the SIS epidemic on networks 

Lloyd P. Sandersp Dirk BrockmannP and Tobias Ambj6rnssoii3 

We provide a closed form perturbative solution to the M-city network SIS model using the trans- 
port rates between cities as a perturbation parameter. We calculate explicitly the first-order per- 
turbation, and quantify that validity of the approximation with respect to the number of initially 
infected and the mobility rates between cities. 

To further corroborate the results, we simulate and analyze a highly infectious SIS outbreak in 
New Zealand (NZ, Aotearoa), where we illustrate how the first-order expansion compares to the 
numerical solution, when we recompose NZ into two subpopulations: Auckland ( Tdmaki Makaurau, 
containing over one third of the total population), and the remainder of NZ. In this comparison we 
use documented air traffic data, and isolate the country from international travel. We demonstrate 
that within the first-order solution's range, it compares well to that of the numerical simulations. 



INTRODUCTION 

In mathematical epidemiology the susceptible- 
infected-susceptible (SIS) model is the one of the 
most elementary compartmental models, receiving 
consistent attention (the specifics of which are given at 
length below) since the seminal work of Kermack and 
McKendrick [1|. Put simply, this model partitions a 
large, homogeneous population into two compartments: 
susceptible and infected. The population is a well mixed 
system, and the birth and death rate are neglected. 
When a susceptible person comes into contact with an 
infected person, there is some finite probability that the 
susceptible person may become infected. An infected 
person will lose their infection after some typical time, 
returning back to the susceptible health state. 

The simplicity of this model and its ability to charac- 
terize the main motifs of viral infections where recovery 
does not assure immunity (for example Gonorrhea [2] or 
Chlamydia [3||), has allowed for extensive research. Due 
to the mathematical tractability of the model, many ex- 
tensions have been applied, to include other important 
dynamical factors; recent investigations include: analyz- 
ing the integrability of the SIS model with so-called vital 
dynamics [4| (explicit inclusion of birth and death rates, 
where new-borns are usually birthed into the suscepti- 
ble compartment); understanding time-delayed infection 
periods [5]; and how intrinsic noise on the transmission 
rate may affect the dynamics [£]. Although these models 
have focused on deterministic mean-field approaches as 
we shall herein, there is also a surge in converting the 
models over to their stochastic counter parts [7], and an- 
alyzing different properties of the system, for a recent 
example see [8|. 

Over the last two decades the literature has seen a pro- 
fusiojL of graph/network theory research (for a review, 
see [9| ) , from which there has been intense focus on how 
this vogue field may contribute to the understanding the 
large scale realism of human mobility [id . Illj . and in 
turn comprehending the etiology of epidemics on these 
systems jl2|. Within approximately the same period 



people have incorporated the concept of metapopulations 
[13| (a population of populations where in each, mean- 
field equations suffice to describe the system dynamics) 
into theoretical epidemiology |14| : where interesting re- 
cent examples include the work by Lund et. al. [15| 
in investigating the effects of heterogeneity on epidemic 
spreading in metapopulations; or the work by Colizza 
and Vespignani on the invasion threshold in heteroge- 
neous metapopulation networks [l6|. 

These similar directions of implicit/explicit spatial 
structure incorporation have naturally been applied to 
SIS models. Of late, network theory has been applied 
to the dynamics of an SIS virus on random networks 
17| , how non-Markovian infection changes these dynam- 
ics on a network [18], or simply exploring the full param- 
eter space of the SIS dynamics on a network via large 
scale numerical simulations |19|. Concerning metapop- 
ulations, Arrigoni and Pugliese considered a stochastic 
Af-city, TV-inhabitant SIS model [20] and Ball analyzed 
the stochastic and deterministic SIS model within a pop- 
ulation of households [2l|. Essentially the current state 
of the research has implied two salient points: com- 
putational power is easily accessible, and the network 
epidemic modeling is not readily amenable to classical 
mathematical tools. In this manuscript we hope to ad- 
dress these topics, whereby we amalgamate the canon- 
ical deterministic SIS model with the current impetus 
toward network/metapopulation modeling through the 
use of perturbation theory, to quantify the effects of hu- 
man mobility on an arbitrary network through analytical 
calculations, substantiated by the equivalent numerical 
simulations. 

Within the following section we review the ana- 
lytics of the canonical single-city SIS model, after 
which we segregate the population into an M-city net- 
work(/metapopulation), whereupon the SIS infection is 
introduced separately to each. We present a closed-form 
recursive perturbative solution to the network model, 
therefrom we calculate explicitly the first-order pertur- 
bation leading to our study's main result Eq. (|20|). Sub- 
sequently we compare our result to a test case scenario 
using real-world population and air traffic data. We then 



discuss the benefits and limitations of the model and 
where this work may be applied and built upon. 

SIS MODEL 

In this section we will describe the equations which 
govern the single- and M-city models and perform ana- 
lytical analysis where applicable. Then, the results are 
compared to numerical simulations. 

Canonical single-city SIS model 

The single-city SIS model considers a large, well mixed, 
population, N, in some closed environment where the 
mortality and birth rates are neglected. The populace is 
divided into two compartments: susceptible, S] and in- 
fected, /; where N = S + I — constant. Susceptibles may 
become infected from contact with the infected at a rate 
/3, and the infected compartment of the population will 
lose constituents at a rate 7. The mean- field dynamics 
of each state is then described by the set of equations: 



dtS{t) = -^SI + ^I, 



dtl{t) 



N 



SI - jl, 



(1) 



(2) 



where initial conditions are: S(t = 0) = 6*0 > 0, and 
I{t = 0) = /q > 0. One immediately notes that 
dt [S{t) + I{t)] — 0, which then ensures that the total 
population, N , is constant for all time. 

Solving Eqs. ^ and ^ analytically, the solution - 
a logistic growth function - is well known (for example 
see J2|), but is included here for completeness. As this 
is a closed system, it is sufficient to consider only the 
solution to the infected population - since S{t) ~ N — 
I{t). Therefore, we can recast Eq. 1^ into: dtl{t) = 
I3{N — 1)1 /N — 7/, and make a transformation: y ^ I^^ 
(such that 7^0), where dy/dl — —P. Therefore Eq. 
^ can be written as dty{t) — ~xy + Pl^^ where x — 
/3 - 7. The solution is thus y{t) = Ce'^* + PI{Nx), 
where C ~ yo — (3/{xN), and yo = I/Iq. Reversing the 
transformation, we find the solution to /(t), 



Iit) = 



1 + Ve-xt -' 



(3) 



where Iqc = X-^//3, is the stable/endemic state of the 
infected population; and V — loo /la — 1- With respect 
to Eq. ([3]), for an epidemic to take place (i.e. some 
finite fraction of the population remains infected in the 
long-time limit), we require the basic reproductive ratio: 
^0 = P/l to satisfy: i?o > 1 (which will be assumed 
henceforth) . 

To incorporate a spatial component to the model let 
us consider an arbitrary network of cities. 



Perturbative solution to the Af-city SIS model 

We here further the result found in the previous sec- 
tion through incorporation of an implicit spatial compo- 
nent, by stratifying the large population into M subpop- 
ulations (caricatured by Fig. [T]). Within each subpop- 
ulation. which may be regarded as a city (but just as 
readily considered to be a community, town, or country), 
the same assumptions stated in the canonical model still 
hold: the population is sufficiently large, and well mixed. 
One then allows for mobility between the nodes (cities) 
on the network, where the fraction of persons traveling 
from city j to i per unit time is given by the transport 
rate, uJi^j, where detailed balance |12l | is assumed: 



Wj^iNi = Wi^jNj, 



(4) 



guaranteeing the total population Ni of city i is constant 
for all time. 




FIG. 1. Caricature of an arbitrary network of communities, 
infected with an SIS type virus, described by Eqs. ([Sjl and ((6|, 
where gray denotes susceptible people, orange, the infected 
and the arrows show possible mobility between nodes. Ni is 
the total population of city i, and (jjj<_i is the travel rate, the 
fraction of the population that travels from city i to j per unit 
time. 



Upon each city, an SIS virus is introduced, where the 
i'^ city has an infectivity rate of f3i , and recovery rate 7^ . 
The mean-field set of equations that then describe the 
dynamics are given by 



N,, 



M 



dtSi{t) = -j^Sih + ^ih + e^ {^t^jSj - Wj^iSi) , 

(5) 



A 



M 



dth{t) = J^^i^^ ~ 7i^j +s'^{'u^t^]I] - Wj^ih) ■ (6) 

i=i 

where we have defined LJi-i-i — 0, and e(= 1) is a conve- 
nient parameter introduced here for counting perturba- 
tion orders (see below). Sunrming Eqs. ^ and ©, we 
find: dtNi = 0, as required by detailed balance. 



To begin the derivation of the perturbative solution, 
we first assume the influx and outflux of citizens from 
a given city is small (defined quantitatively later), from 
there we can define a perturbative solution in terms of 
the travel rates, namely 

oo oo 

hit) ^Y.''ll'\t) = /f )(i) +^e'^7W(i), (7) 



fc=0 



fc=i 



where I^ (t) is given by Eq. ([3]) , with replacements (3 -^ 
f3i and 7 — ?> 7i (and therefore x ~^ Xi)- Explicitly: 



r(0) 



irit) = j 



/™ 



+ y,e-x.* • 



(8) 



Similarly we can define the perturbative solution to the 
number of susceptibles in city i a.s Si = X]fc=o^'^'^i ' 
which then implies, due to detailed balance, that 



r(fe) _ 



'S. 



(k) 



for fc > 1. 



(9) 



Substituting Eqs. dH) and ^ into Eq. dH) we find 



d.Y^e'^I, 



kAk) 



A 

N, 



E 

k.l 



e^e^ 



r(fe) 



:E^'=^'^+E -»-.E-'"^^ 



-fc+1 A^) 



E 



fc+l r(*:) 



Equating factors of e , we find for fc = 0, that 

2 



dj^'' 



Xtl. 



(0) 






(0) 



(10) 



(11) 



which is equivalent to Eq. ^, and whose solution 
is therefore given by Eq. ([5]), with appropriate i- 
dependencies of the virus parameters. For fc > 1 we 
obtain our formal perturbation equations 




2(3d, 



(0)- 



N,, 



T 



(fc) 



r(*:-l) r 




(k-k') 



{k-1) 



(12) 



Thus, we have formally converted the non-linear prob- 
lem in Eqs. ([5]) and ([H]), into a set of inhomogeneous, 
linear equations - Eq. p2p - with time dependent coef- 
ficients. The time dependence of these coefficients en- 
ters only through the known quantity, /°(i), whereas the 
right-hand side depends recursively on the previous per- 
turbation orders. 

We proceed by expressing a formal solution to 
Eq. p2p though employment of the integrat- 

ing factor method. Firstly, we define the so- 

called integration factor: exp(i?i(t)), where Bi{t) = 



r(0)^ 



Lth 



dt' . Then the formal solu- 
tion to the fc"^-order perturbative term is /| ' [t) = 
exp(-Bi(t)) J^exp{B,{t'))g,{t')dt' + G , where from 

the initial conditions: /| ' {t = 0) =0, we have G — 0. 
The function gi{t) is defined as 

.(fc-l) ,,_ r(fc-l)' 



9i 



\t) 



N^ 



+ E 



wj^j/j' 



- UJ 



j^^h 



(13) 



Interestingly, we are able to calculate Bi{t) explic- 
itly. Using Eq. ([5]), we can write Bi{t) = —Xit + 
{2l3Jooa)/{Nr)J^{l + V^e-^^*')-^dt'. We solve this to 
yield the solution 



B^{t) = In 



^Xit 






(14) 



Using the solution for Bi{t) and Eq. ([T^ and substi- 
tuting this into the formal solution given, we explicitly 
obtain the fc'^ order perturbation, namely 



4'\t) 






(15) 



With this closed form expression, we are able to calculate 
any order perturbation we require, recursively. Namely, 
starting from the zeroth-order solution, Eq. ^, we can 
insert this into Eq. (jlSp , the result of which is then input 
into Eq. ([TS]) to find the first-order perturbation (shown 
explicitly in the following section) . If one is to find the 
next order, one merely uses the first-order result in place 
of the zeroth-order solution to, following the outlined al- 
gorithm, arrive at the second-order perturbation. This 
operation may be repeated until the desired number of 
orders are achieved. Then the orders are summed, viz. 
Eq. ([7]), to gain the final solution to the infected popu- 
lation contained in city i. 

Explicit first-order perturbation 

Let us now calculate the first order perturbation term, 
for which fc = 1. Then function g\ , see Eq. (fT3|) . is ex- 



r(0) 



W,-^,/' 



(0) 



such that 



plicitly: gf\t) ^J2j[^^^3 1 J - ^j^^^i y> 
Eq. (fT5| . using Eq. ([8]), becomes 

-X^t 

y , [i^i^jQij - ^j^iQu] (16) 



j(i) ^ e 



(1 + Vie- 



Xit\ 



where 



Q 



ij — -tooj 



/oo., / e^'* ^^ ., . / dt'. (1 



1 + VjC-xA' 



(17) 



The quantity Qij may be expressed in terms of hyperge- 
ometric functions |22| . 

For the scope of this manuscript, let us analyze the 
case where we shall assume that all virus rate parameters 
are independent of the city, namely j3i = /3j = /3 and 
7i — 7j- — 7. Explicitly evaluating Eq. ([17]), we find that 



'ooj 



{V^-Vjf ^^(l + V,e-^' 



V, V l + ^j 



^Xt{2V, -~Vj) + l- e^' 



(18) 



Due to detailed balance, Eq. (|3]), we know that 
using this relation, we can write 



LO, 



-3 -'-00,3 



-i-'oo,? ; 



out the first order perturbation, Eq. (|T6|) . explicitly as 



r(i) 






^{V^~V,)X 



xt 



V^ ~ V3 



In 



1 + Vj-e-x* 

~TTv~ 



(19) 



Summing this with the zeroth-order solution we reap the 
first-order perturbative approximation to the number of 
infected in city i at time t: 



Ht) = Y 



+ y,e-x* 



-xt 



(1 + V^e 



-xt 



■E^ 



''-{V^-V,) 



xt- 



V, 



Vi 



-^m 



1 + Vj-e-x* 
1 + K- 



0{w}^,] 



(20) 



r 



where 



and 



V^ 



lo. 



^00.7 






The first-order approximation to Ii (t) , namely. Eq. (P0)l , 
is valid when the perturbation due to mobility is "small". 
To quantify explicitly the validity of Eq. ([20]) we intro- 
duce the validity criterion: 



C, 



X 



-{V^~V,) 



lo.i 



lo.. 



(21) 



For a given system, when ^ ■ Cj^i « the zeroth-order 
solution is valid for Ii{t); when ^ ■ Cj^i > 1 the first- 
order solution breaks down. In particular we point out 
that, in fact, besides the travel rate ujj<^i (in units of /3), 
also the fraction of initially infecteds enter Eq. (|2ip in a 
non-trivial way |23| . 



SIS Epidemic in New Zealand {Aotearoa) 

We proceed by measuring the performance of our an- 
alytical expression, Eq. (|20p . with realistic population 
and air transport data through a test case scenario, an 
SIS epidemic in New Zealand (NZ). Firstly, we suppose a 
highly infectious virus (/3 = 0.25 day" , 7 = 0.025 day" , 
Rq = 10) has established itself within NZ, and assume 
the virus is only spread city to city via the air traffic 
transport network where there is no transport interna- 
tionally with NZ (which one could think of as a quaran- 
tine measure). Then, to understand how the country's 



largest city, Auckland {Tamaki Makaurau), iVAuck, is af- 
fected/or affects the remainder of NZ, A^rem, we recom- 
pose NZ into these two subpopulations (thereby assuming 
that each of these populations are large - whose explicit 
populations are stated in Fig. [5]- and well mixed). Real- 
istic transport rates are Wiom-i-Auck = 6.0 x 10""^ day" , 
and WAuck^rcm = 3.1 X 10"'^ day" [24]. We set ini- 
tially the infected to /o,Auck = 0.01 x iVAuck, smd /o,icm — 
0.026 x A^rem, for Auckland and the remainder of NZ, re- 
spectively. With these parameters we note that validity 
criterion, Eq. (PT|) . for Auckland, is Crem<-Auck = 1-45 
and for the remainder of NZ is CAuck<-rem = 0.74. 

We solved the dynamics of the individual infected pop- 
ulations of this system numerically (by the Runge-Kutta 
^th order algorithm), given the parameters above and 
compared the numerical results to our result in Eq. (I^Ul) . 
The results are plotted in Fig. [2l and are compared with 
the scenario that all airports are shut down, which is 
also a measure of how effective the the zeroth-order solu- 
tion, Eq. ([5]), is at estimating the time-evolution of the 
epidemic. We note that our analytical result conforms 
well to the numerical result, where the absolute residuals 
(= l^numcr.(i) " Ii(t)\, where the former is the numerical 
result, and the latter our result from Eq. (HO])) between 
both are given in the inset. The zeroth-order, "no travel," 
solution performs poorly at estimating the interim dy- 
namics of the epidemic, especially underestimating the 
fraction of infected in Auckland (at its peak residual, 
this result is off by approximately 0.03 x A^Auck = 44580 
people.). 

To assess the scope of our perturbation result in this 
system, we explore the phase space of the fraction of ini- 
tially infected for each subpopulation - Fig. [3] - with 
respect to Auckland (the corresponding phase-space plot 
for the remainder of NZ is qualitatively comparable.). 
We linearly vary the fraction of infected from 10"'' -^ 



2 X 10~^, and for each pair we calculate the perturbation 
approximation, and the numerical solution. The abso- 
lute residuals (see Fig. [51 inset) between the numerical 
and perturbation solutions are calculated; the area under 
which is found. This area is divided by the area under 
the numerical solution to give the absolute residue area 
ratio. These ratios are plotted as a "heat map" in Fig. 
[HI where the lighter the shade, the closer our analytical 
result is to the numerical solution. When the fraction 
of initially infected for each subpopulation, compared to 
the other, is near equal the absolute residue area ratio 
is relatively small, and therefore a good approximation. 
This is predicted by our quantitative benchmark, the va- 
lidity criterion, whose contour lines have been overlaid 
upon the heat map. The usefulness of the first-order ap- 
proximation lies in the off-diagonal circumstances, and 
for one such instance is shown, via a white circle, on Fig. 



[31 which is the pair of the initial fractions used for Fig. [31 
These initial conditions demonstrate the outer limits of 
the approximation, where the zeroth-order has effectively 
failed, and the first-order is able still to approximate the 
numerical solution. 

If we instead push the approximation past the recom- 
mended parameter sets, we can observe its breakdown. 
We present this in Fig. [H using the initially infected 

values: (/o.Auck,/o,rem) = (0.0025 X iVAuck,0.05 x iV,.cm)- 

This point is shown on the phase-space plot. Fig. [31 as a 
white square. 



DISCUSSION AND CONCLUSION 

Within this manuscript, we have derived a closed form, 
recursive, perturbative solution to an SIS epidemic on an 
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FIG. 2. The infected fraction of the populace over time 
(due to an SIS epidemic) for each city for an internationally 
quarantined New Zealand, apportioned into two subpopula- 
tions: Auckland, and the remainder of the country. Indi- 
viduals (both susceptible and infected) may migrate between 
these two subpopulations only via commercial air transport; 
and through documented air transport data [25||, we find the 
transport rates are effectively small - viz. Eq. (|21|l . One 
clearly notes that in both populations, the first-order per- 
turbation, Eq. (|20p (solid red/black line), approximates the 
numerical solution well (solid gray line), compared to the 
zeroth-order solution, Eq. (|H| (dashed red/black line), the 
case where infected migration is not considered. Transport 
rates (see main text) are: tJrem<-Auck = 6 x 10^^ day~^, 
and ojAuck^-rem = 3 x 10"'' day~^ [2J|. Validity criteria, 
Eq. (|2H) . for each subpopulation are: (7rem^Auck = 1-45 and 
CAuck*-rem = 0.74 (scc also Fig. 13). INSET: Absolute residu- 
als of the first-order perturbation and the zeroth-order solu- 
tion, with respect to the numerical solution (defined explicitly 
in the main text). Auxiliary parameters: population sizes (for 
the same year as the air trafiic data), A'Auck =1486000, and 
Afi-om =2919300 [26|; virus parameters are: /3 = 0.25 day~^, 
7 = 0.025 day~^. Initially infected populations for the two 
cities are: /q Auck = 0.01 x iVAuck, and Jo, rem = 0.026 x A'^rem- 



FIG. 3. Exploration of the phase space for the fraction of 
initially infected in the quarantined NZ system, with respect 
to Auckland. The lighter the shade, the closer the pertur- 
bative solution for the infected population in Auckland is to 
the numerical solution (explained below). Linear alteration 
of the fraction of initially infected for the Auckland popula- 
tion is shown on the ordinate, and for the remainder of NZ 
population on the abscissa (all other system parameters are 
kept constant - see Fig. [2]). For each set of initial popula- 
tions, we compare our perturbation result, Eq. (|20p . to the 
numerical solution. The absolute residues for Auckland (see 
Fig. [2] inset) between these two results are calculated, then 
the area under which is measured and divided by the integra- 
tion of the numerical solution for the same time-frame. The 
absolute residue area ratio is plotted as a heat map, where 
the shade (see color bar above) denotes the magnitude of the 
value. The initial conditions used for Fig. [2] are marked with 
a white circle in the lower left quadrant. An instance when 
the first-order perturbation breaks down is highlight via a 
white square, and whose expHcit temporal evolution is shown 
in Fig. [4] Superimposed on the initial condition heat map is 
the corresponding contour plot of the validity criterion, Eq. 
(|2ip (where C denotes Crem<-Auck)- The phase space heat 
map for the remainder of NZ is qualitatively similar. 



arbitrary network, stated in Eqs. ([T^ and (|15p . We have 
proceeded to explicitly calculate the first-order perturba- 
tion to the population of infected persons in the z*'' city 
as a function of time, to wit, Eq. ^IU\\ . and then pro- 
vided a quantitative benchmark under what conditions 
this solution is accurate: the validity criterion, Eq. (PT|) . 

To verify our derived results we simulated an SIS epi- 
demic on an internationally quarantined New Zealand. 
This comparison served a two-fold objective; firstly the 
use of documented air traffic data [25| showed that in 
this medium of transport, our base assumption: trans- 
port rates between communities are small, is indeed rea- 
sonable and accurate as a perturbation parameter (see 
main text for further discussion, caveats, and explicit 
numbers). Secondly, when the system parameters give 
a validity criterion near or less than unity, the first-order 
solution, Eq. (|20p. approximates the numerical solution 
well, as shown Fig. [2l The range of usefulness of the 
validity criterion and its dependence on the fraction of 
initially infected is explored in Fig. [3l and a set of pa- 
rameters for when the approximation breaks down is pro- 
vided in Fig. m 

The derivation and the subsequent verification through 
the numerical simulations assumed detailed balance and 
has leaned on the assumption that a finite fraction of 
each city is initially infected; the derivation then seeks 
to describe temporal evolution of the epidemic. In this 
case, this would apply to the reasonable scenario where 
the authorities of a region would look to understand the 
dynamics of the epidemic once their community of cities 



Numerical solution — 

Auck. (pert.) — 

NZ(pert.) — 

Auck. (no travel) — 

NZ (no travel) -^ 



a 
o 



-a 



1 

0.9 
0.8 - 
0.7 
0.6 
0.5 - 
0.4 
0.3 
0.2 - 
0.1 




FIG. 4. Demonstration of the breakdown of Eq. (|20|) . our 
first-order perturbation, with initial infection parameters cor- 
responding to a large validity criterion. In this system the 
initially infected for Auckland and the remainder of NZ are: 
/o.Auck = 0.0025 X A'^Auck and Jo, rem = 0.05 x A''rem respec- 
tively (this pair is marked as a white square in the phase- 
space for Fig. [31). The corresponding validity criteria are: 
Crem^Auck = 9.14, CAuck^rem = 4.65. For the remaining sys- 
tem parameters and explanation of inset see Fig. [5] caption. 




is found to be infected. This of course shifts the prob- 
lem from the usual discourse found within the literature, 
where the infection is present initially in one population 
and spreads, such as [27|. It would be of interest to un- 
derstand how higher order perturbation terms may reg- 
ularize our first-order approximation. 

One of the striking benefits of the solution is that one 
has an analytical expression for Ii(t) at all times and as 
such naturally out-performs usual investigative methods 
of numerical integration. In this way this solution can be 
used to gauge parameter sets of large numerical simula- 
tions. 

The derivation makes no assumptions on the type of 
the network, whether it be a real-world network, regular 
lattice, a random E-R graph, or a scale- free network [9|. 
It has also only assumed that the population of each node 
is large enough such that the mean- field nature of Eq. ^ 
is true. Therefore the nodes may be seen as communities 
or countries, rather than only cities. This generality is 
advantageous for future investigations of SIS epidemics 
on complex networks/metapopulations as this work may 
be used in parallel as a confidence measure. 

These calculations were built upon the assumption of 
large populations, where mean-field approximations are 
valid. A natural extension would be the effect of stochas- 
ticity for low populations. This may be found through 
an analytic perturbative solution of the associated mas- 
ter equation akin to that defined for an SIR epidemic in 
the work of Hufnagel et al. 27|. 

Although the analytics have been developed under the 
guise of epidemic modeling, this mathematical frame- 
work may be conveniently adopted by other interdisci- 
plinary fields with population growth and metapopula- 
tion structure, for example Theoretical Ecology and the 
concept phenomenon of island colonization [28| . 

In conclusion, we hope the mathematical framework 
determined herein will shift part of the academic inter- 
est of epidemics on networks from large scale numerical 
simulations back to the bedrock of analytical analysis. 

The authors would like to thank Bo Soderberg, Erik 
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son for stimulating discussions regarding the work herein. 
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